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In this paper we discuss a new method which can be used to obtain arbitrarily accurate analytical 
expressions for the deflection angle of light propagating in a given metric. Our method works by 
mapping the integral into a rapidly convergent series and provides extremely accurate approxima- 
tions already to first order. We have derived a general first order formula for a generic spherically 
symmetric static metric tensor and we have tested it in four different cases. 
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I. INTRODUCTION 



According to general relativity the trajectory of a ray of light which passes close to a mass distribution departs from 
, being a straight line. The amount of the deflection of the light depends upon the mass and can be quite large for light 
7-H ■ passing very close to a massive compact body, such as a black hole. The study of gravitational lensing under such 
[ conditions, also known as "strong gravitational lensing", has received wide attention in the recent past: for example 
strong gravitational lensing in a Schwarzschid black hole has been considered by Frittelli, Kling and Newman Q 
' and by Virbhadra and Ellis 0; Virbhadra and Ellis 01 have later treated the strong gravitational lensing by naked 
singularities; Eiroa, Romero and Torres Q have described Rcissncr- Nordstrom black hole lensing, while Bhadra has 
considered the gravitational lensing due to the GMGHS charged black hole ; Bozza has studied the quasicquatorial 
gravitational lensing by a spinning black hole ; Whisker and Eiroa @| have considered strong gravitational lensing 
PsJ ' by a braneworld black hole; still Eiroa |3| has recently considered the gravitational lensing by an Einstein-Born-Infeld 
' black hole; Sarkar and Bhadra have studied the strong gravitational lensing in the Brans-Dicke theory finally 
Perlick has obtained an exact gravitational lens equation in a spherically symmetric and static spacetime and 
' used to study lensing by a Barriola-Vilenkin monopole and by an Ellis wormhole. 
O . Different strategies have been used to evaluate the effects of strong gravitational lensing: for example, Bozza [l^ 
O^' has introduced an analytical method which allows to discriminate among different types of black holes: the method 
, is based on a careful description of the logarithmic divergence of the deflection angle (the photon sphere); Mutka and 
Mahonen 0, and Belorobodov have derived improved formulas for the deflection angle in a Schwarzschild 
^ [ metric; more recently, Keeton and Fetters [l6j have also developed a formalism for computing corrections to lensing 
observables in a static and spherically symmetric metric beyond the weak deflection limit. 

The purpose of this paper is to present a new method which can be used to calculate analytically and systematically 
the deflection angle in a static spherically symmetric metric. Originally this method was devised by Amore and 
collaborators |17l to obtain analytical formulas for the period of a classical oscillators: our method works by 
converting the integral which needs to be calculated into a scries dep ending upon a variational parameter. Such 
procedure is inspired by the Linear Delta Expansion (LDE) method and by Variational Perturbation Theory [23| • 
For certain values of the variational parameter, the series obtained is proved to converge to the exact result, while at 
finite orders a particular value of the parameter can be chosen using the Principle of Minimal Sensitivity (PMS) 
to minimize the error. Fully analytical results, which do not correspond to a perturbative expansion in some small 
parameter, are obtained. 

The paper is organized as follows: in Section^ we describe the method and obtain a general first order formula, 
which is valid for an arbitrary metric tensor; in Section [llll we discuss the convergence of the method and provide an 
estimate of the rate of convergence, which is proved to be exponential; in Section IIVI we apply our formula to four 
different metric tensors and discuss the precision of our approximation, comparing it with the available results in the 
literature; finally in Section^we draw our conclusions. 
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II. THE METHOD 

We are interested in the general static and spherically symmetric metric which corresponds to the line element 

ds^ = B{r)dt^ - A{r)dr^ - D{r)r^ {d0^ + sin^ 8 d(f>^) (1) 

and which contains the Schwarzschild metric as a special case. We also assume that the flat spacetime is recovered 
at infinity, i.e. that limr_»oo fir) = 1, where /(r) = {A{r), B{r), D{r)). 

The angle of deflection of light propagating in this metric can be expressed by means of the integral |22| 



A0 = 2 / ^Air)/Dir) 



\ 



rnj D{ro) B{r) 



--^ , (2) 

r 



where ro is the distance of closest approach of the light to the center of the gravitational attraction. 
If we perform a change of variable, z = tq/t, and define the function 



^ ' - A{ro/z) A{ro/z)B{ro/z)Diro) ^ D{ro) 



we can write eq. (0) in the form 



= 2 /' , - ^ . (4) 

io ^V{l)-V{z) 



Notice that t = y/2 dz/ y/V{l) — V{z) is the time spent by a classical oscillator moving in a potential V{z) for 
passing from z = to the inversion point located at z = 1; in a flat spacetime, where A(r) ~ B{r) ~ D{r) = 1, 
V{z) reduces to the familiar harmonic oscillator potential and the deflection angle identically vanishes. The integral 
of eq. @ can also be performed exactly in the case of the Schwarzschild metric and of the Reissner-Nordstrom 
metric 0. In both cases the exact result is expressed in terms of elliptic integrals. However, in more general cases 
the integration of eq. Q cannot be done exactly and one typically resorts to an expansion around the flat metric (of 
course such expansion can also be used in cases where exact results are available, to avoid dealing with complicated 
special functions). In the case of the Schwarzschild metric, for example, this approach yields a perturbative series in 
powers of GM/rQ, whose leading term A<j>o = 4:GM/ro was first obtained by Einstein. A big disadvantage of this 
approach is that the validity of the expressions obtained in this way is restricted to large distance/ weak field regime: 
for example, the exact solution to lensing in the Schwarzschild metric provided by Darwin possesses a singularity at 
ro = 3GM, known as photon sphere, and is clearly out of reach in a perturbative approach. 

We will therefore pursue a different approach to deal with eq. I0J, which is not based on a perturbative expansion 
and which we will prove to be capable to describe very accurately the physics of our problem. The method that we 
propose has been devised by Amore and collaborators 0, 0| to obtain precise analytical formulas for the period of 
classical oscillators. In a recent work, the method has also been used to obtain analytical expressions for the spectrum 
of quark-antiquark potentials |24j . A similar technique has also been applied by Amore to accelerate the convergence 
of certain series (such as the Riemann and Epstein zeta functions) j25il , which can be used in the calculation of loop 
integrals in finite temperature problems occurring in field theory [26|. 

We now briefly describe our method. In the spirit of the Linear Delta Expansion we interpolate the full potential 
V{z) with a solvable potential Vo(z)^: 

Vs{z) = Vo{z) + 5{Viz)-Vo{z)) . 

Depending on the value of S one will obtain the original potential ((5 = 1) or the solvable potential {S = 0). In 
general Vo{z) will depend upon one or more arbitrary parameters, which we will call A; in fact we will assume the 
form Vo{z) = Xz^ . With this definition wc write the deficction angle as 

. , dz 

Jo VM^) - Voiz) + 5{Vil) - Viz) ~ Voil) + Vo{z)) 



^ By solvable we mean that integrals z"dz/^/V{l) - V{z) can be performed analytically. 
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which clearly reduces to the original expression for S = 1. 
After introducing 



A(z) 



Eo-Vo{z) 



one can write (O as 



A(j)s 



dz 



^E„-Vo{z) Vl + <5A(2) 



(6) 



(7) 



Provided that |A(z)| < 1 for < z <1 one can expand eq. (UJ in powers of i5 and obtain a series (after performing 
the integrals) which converges to the exact result. As discussed in this condition requires that A be greater 

than a critical value, A > Ac: in this case one obtains a family of series which depend upon A and which all converge 
to the exact result, which is independent of A. However, if the series is truncated to a finite order, the partial sum 
displays an artificial dependence on A: such dependence can be minimized by applying the Principle of Minimal 
Sensitivity (PMS)[23: 







(8) 



having called Ac/)'^^ the partial sum to order TV. 

Notice that the solution to this equation selects the value of A where the series is less sensitive to changes in A itself: 
this value of A selects the series with the optimal convergence. In a large class of oscillators was studied using 

this method and it was found that our PMS series has an exponential rate of convergence. 

However, since it was observed that the first order results are quite precise, we focus our present effort in obtaining a 
first order formula, which is valid for a generic spherically symmetric static metric tensor corresponding to a potential 



V{z) 



CO 

n=l 



2 



After expanding to first order we obtain 

A0(i) . 
where 

A{z) 

The deflection angle can now be written as 

A0(i 



dz 



T 

n=l 



fc=0 



+ z 



2^/A 



A3/2 Z-^ 

n—1 



n-1 

E 

fc=0 



where we have defined 

h = 



^'^ sin^ 



-de = 



r((fc + i)/2) 



r(fc/2 + 1) - 



[r((fc + i)/2)f 

T{k/2) 



(9) 



(10) 



(11) 



(12) 



(13) 



Notice that eq. l(T^ can be expressed in terms of the "transformed" potential 



p{z) = ^ w„ri„z" 



(14) 



where 



n-l 

E- 

A:=0 



-r(n/2 + 1/2) 
r(n/2) 



(15) 
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As expected our first order result depends upon A and we must use the PMS to obtain the optimal value of A: 

^PMS - ^ ' I -"-"J 

Once this value is substituted inside eq. p2|) we obtain our first order result 



2p{l) 



1 



(17) 



which in our opinion is the most valuable formula contained in this paper. Notice that because of the form of eq. p7|l 
our approximation does not correspond to a perturbative expansion in some small parameter, as it will be clear in 
the next section. 

III. CONVERGENCE 

In this section we discuss the convergence of our method: we wish to prove that our procedure provides series with 
an exponential rate of convergence. As we have mentioned in the previous section we can write the expression for the 
deflection angle in a power series in A(z) as 



^^ = 2f^J^= v(?!i4l!lA"(.)-.^VA„, (18) 



if the condition |A(2)| < 1 is fulfilled in the region of integration, z E (0, 1). 

Let us now call A^ax the maximum value of |A(z)| in the region of integration; we can therefore write 

|A„| < A-^a. ■ (19) 

For large values of n one can approximate the factorial and double factorial in this expression with the corresponding 
asymptotic series, (2n — l)!!/2"7z! sa l/^/nn and therefore obtain 

|A„| < , (20) 



which confirms that the series converges geometrically. Notice that since A^ax is defined in terms of the original 
potential V{z) and of the interpolating potential Vo{z), it will depend upon the arbitrary parameter A: we therefore 
expect that the optimal value of A, obtained at a finite order using the PMS, will be such that at large orders A^ax 
assumes the smallest value possible. Alternatively, one could think of lowering the value of Amax by choosing a 
potential different from the simple harmonic potential, Vo{z) = Xz^, discussed in the previous section: indeed the only 
limitation that we have provided over Vo{z) is that it is such that the integrals contained in the scries for A(j) can be 
performed analytically. This strategy, although possible, is not followed here because the increased complication in 
the form of Vq{z) would necessarily reflect in a complication of the formulas obtained and in the drawback of obtaining 
approximations in terms of special functions. Clearly, such a procedure should also be investigated in future works. 
We will now provide an estimate of Amax' under the assumption that A{z) is a monotonous function, we have that 

A„,, =max{|A(0)|,|A(l)|} . (21) 

The convergence of the series requires that A„iax < 1, which is fulfilled for 



A > Ac = niax-^ -^nt;„,-^t;„ ^ . (22) 

I ^ n=l ^ n=l J 

On the other hand it is easy to convince oneself that the minimal value of Amax is obtained when the condition 
A(0) = -A(l) is fulfilled, i.e. when 



71=1 



^ 00 
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corresponding to 



Er=it(i + «/2) 



The reader should also notice that our scries cannot be used when \max < Ac, because the condition l^max < 1 
cannot be obeyed. 

We can easily test our results over the Dufhng potential V{z) = 2^ + z'^; in this case we have that Xmax and Xpi\is 
coincide {XpMS = 5/2). Corresponding to this value we have Amax = 1/5 and we obtain a rate of convergence which 
is stronger than r„ w (l/5)"/-\/7rn, is in good agreement with the rate observed fitting the behavior of the series up 
to order 10, r„ w 0.028 x 0.17". 



IV. APPLICATIONS 



We consider in this section four applications of the formula H17|l obtained in the previous section. In the first 
two cases explicit formulas for the exact results are known due to Darwinp^ and Eiroa and collaborators 0; in the 
last two cases we consider the metric of Janis-Newman-Winicour and the metric of a charged black hole coupled to 
Born-Infeld electrodynamics, for which no explicit formula is available. 



A. Schwarzschild metric 



B{r)=A-'{r)=[l-^^] , D{r) ^ 1 . (25) 



Our first application is to the Schwarzschild metric, which corresponds to 

'-) 

Here M is the Schwarzschild mass. The angle of deflection of a ray of light reaching a minimal distance rp from the 
black hole can be obtained using eq. Q . The exact result can be expressed in terms of incomplete elliptic integrals 
of the first kind 1231 and reads 



f(^^,k) -F{^,k)] , (26) 



Acf) = 

where Tq = rO/GM and 

l-r. - 9 I r 9 ^ T - r n " 

(27) 



Iro 


- 2 




+ 6 



6 + T - r 



Although eq. 1)26(1 is exact, it is often valuable to obtain approximations which do not involve special functions. 
Here we will compare our first order approximation, corresponding to using eq. I|17|) , with other approximations which 
have been derived in the literature. 

For example, Mutka and Mahonen 0,^3 have obtained the approximate formula 

Acj)MM = ^ , (28) 



where b = ro\/D(ro)/B(ro) — ro/ \/l — 2GM/ro is the impact parameter. This formula is a natural extension of the 
Einstein formula 

= ^ . (29) 


Beloborodov|l5j has obtained another approximate formula which reads 

4GAf , , 

ro — 2GM 



6 




ro 



FIG. 1; Deflection angle as a function of ro assuming GM — 1. (color online) 



Finally, Keeton and Fetters have devised a systematic approach to deal with integrals as and obtained the 
formula 



A(l) 



KP 



, fGM\ , [GmV' , [GmV , fGM\' 



o 



(31) 



where the numerical values of the coefficients Ai are given in eq.(25) of [iq . 

Using the general equation for the deflection angle to first order, eq. (|17|l . we have obtained the formula: 



corresponding to Ap|jg = ^1 — 8GAf /Trro. 



1 



^1 - 8GM/nra 



(32) 



Despite its simplicity, we can appreciate from Fig. ^ and [3 that eq. (|32|l provides the best approximation to the 
deflection angle, even in proximity of the photon sphere (the singularity): indeed our formula predicts the location 

of the singularity at = 8GM/t: sa 2.55GA'jf, slightly below the exact value ri"^^^ = 3GA/. While the expression of 
Beloborodov puts the singularity at a smaller value of tq, the remaining approximations cither put it in the unphysical 
region (ro < 0) ( Mutka and Mahonen) or fail to produce a singularity (Keeton and Fetters). 

In Fig. 121 we have also plotted the analytical approximation of Bozza which correctly describes the photon 
sphere: our first order formula provides better approximations already for rg > AGM . 

Remarkably our expression works very well also in the opposite regime, corresponding to ro ^ oo; our eq. (|32|l can 
be expanded for ro 3> 1 to give 



AGM 7.63944G2Af2 IQ.2\1AG^AP 



O 



{GM/roY 



which compares quite favorably with the exact asymptotic behaviour of the Darwin solution: 

4GA/ 7.78097G2Af2 i7.i047G^A/3 



ro rj 

In Fig. 13 we also show the magnification (see 



O 



{GM/roY 



M = 



1 



-A0 



db dro J 



(33) 



(34) 



(35) 
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FIG. 2: Deflection angle as a function of ro assuming GM = 1 close to the photon sphere, (color online) 
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FIG. 3: Magnification as a function of the impact parameter assuming a — 1000 and GM — 1. (color online) 



as a function of the impact parameter b. a is the distance between the lens and the source. Once again our simple 
formula provides a very accurate approximation to the exact result over a wide range of values. 

Given the success of our approach to first order, a higher order calculation is not essential, although it is not 
technically difficult^. For example, it is straightforward to obtain the second order formula 



TT (l5^2 _ 1)2 _ i67r (A2 - 1)' - 32 (8A^ - SA^ - 6A^ + 3) 

256A5 



where 



A 



^PMS 



1 - 



8GM 
Trro 



(36) 



(37) 



^ In llTl . for example, the method was applied up to order 100 to calculate analytically the period of an anharmonic oscillator with a 
precision of about 10"®". 
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This formula approximates the deflection angle to a 1% level up to tq w 3.5 GAf , i.e. quite close to the singularity, 
compared to « 7.4 GM of the first order result. 

Using the results obtained in Section Hill we can also estimate the rate of convergence of our series. In this case the 
potential is 

V{z) = z^- (38) 

and 

GM/ro 



2 - 5GM/ro ' ^^^^ 

It is interesting to notice that the condition of applicability of our series. \max > Ac, can be fulfilled only for 
ro > 3GA/, which is the exact location of the photon sphere for the Schwarzschild metric: in other words, our series 
can also describe strong gravitational lensing close to the photon sphere. 



B. Reissner-Nordstrom metric 



The Reissner-Nordstrom (RN) metric describes a black hole with charge and corresponds to 

2GM 



B{r) ^A-^{r) = 



1 



D{r) = 1 



(40) 



As for the Schwarzschild metric the angle of deflection of a ray of light reaching a minimal distance tq from the 
black hole can be obtained from eq. (|2J). Eiroa, Romero and Torres Q have been able to express the deflection angle 
in terms of elliptic integrals of the first kind (see eqn. (A3) of 0). 

It is straightforward to use our general formula to obtain the transformed potential for the RN metric 

^-^^ (41) 



TT AGM 



and thus the deflection angle 



PMS 



SGM 



(42) 



We can compare our formula both with the exact analytical result of Eiroa et al. and with the expressions (47) 
and (53) of [I^, which provide a systematic expansion of the deflection angle in terms of GM/b. In Fig. 01 we have 
plotted the exact solution of Q together with our first order formula and with the expression of Keeton and Fetters, 
assuming G = M = 1 and Q = 1/2^: the reader can appreciate that our simple formula is very accurate even in 
proximity of the photon sphere. 

Notice that our expression reproduces well also the asymptotic behaviour of the deflection angle. In fact, eq. (19) 
of provides the leading asymptotic behavior of Acj), valid for r oo: 



A(j) 



AGM AG^AP 



AGM 



7.78 



G^M^ 



-2.36 -V + 



which can be compared with the asymptotic behaviour of our formula 



AGM 2AG^M'^ 



3tt Q2 



ro 
AGM 



O 



'0 



-O 
1 

ro 
1 

ro 



(43) 



7.64 



G2Af2 



ro 



G^ 

2.36 -V + O 

rl 



(44) 



^ Notice the different definition of Q in Il6l . 
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FIG. 4; Deflection angle for the Reissner-Nordstrom metric as a function of ro assuming GM = 1 and Q = 1/2. (color online) 



Once again we can refer to the results of Section HIH to estimate the rate of convergence of our series. In this case 
the potential is 



V{z) 



2GM 3 ^ 0^ 



and 



(45) 



(46) 



In this case, the condition of applicability of our series. Xmax > Ac, can be fulfilled only for j'q > 



3GA//2 1^1 + i/l — 8/9((5/GA/)^j , which is the exact location of the photon sphere for the Rcissner-Nordstrdm metric 
(see eq.(8) of 0). 



C. Janis-Newman-Winicour metric 



We consider now the spherically symmetric metric solution to the Einstein massless scalar equations [22 
A{r)^{l-b/ry , B{r) = (l-b/r)" , D{r) ^ {1 - b/r)^~'' , 
which reduces to the Schwarzschid metric for v ~ I and for b = GM . In this case we obtain the potential 



which can be expanded around z = to give 
V{z 

1(1, - lH2,y - 1) 
3 V^o 



z'll 



bz 
ro 



1 



_2(i.-l)(l- A)2--iAz- 
?'o ro 



1 - 




-l)(2i. 


-1) 






b 




b ' 












z^ 




ro 




ro 





b_ 

ro 



b_ 

ro 

2v-\ 



Using this expansion inside our formula, cq. H17|) . wc obtain the deflection angle 



(i.-l)(l-;)^'-C((2t^-l);(8i/C+3x ) + 12)-12(;-l)C 



37r(5-l) 



(47) 



(48) 



(49) 



(50) 
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FIG. 5: Ratio between the approximate first order formula and tiie exact (numerical) result in the JNW metric for different 
values of v and assuming h —\. The potential is expanded to order . (color online) 



where ^ = fc/?'o- 

In fig. |Slwe have plotted the ratio A0p|y^g/A(/)ea;act for three values of v, up to small values of tq. Notice that the 
exact result is calculated numerically and that the ratio is close to 1 up to very small values of tq. 
The expansion of eq. lf3n|l around ro reads 



. ^ + (2(6-4. + ^-)^--3(-4 + .)..+ (-4 + .).)g ^ ^ ^^^^3^ 



ro ' 27r '/Q 

and provides a deviation from Einstein's leading order term, Acf) « 2b/ro for ly ^ 1. 



D. Einstein-Born-Infeld black holes 



As a final example of application of our method we consider the propagation of light in a charged black hole coupled 
to Born-Infeld electrodynamics. This problem has been recently considered by Eiroa in [H and corresponds to the 
effective metric 



Mr) = ^^^ , B{r)^^^^j{r) , D{r) 



1 



where 



w(r) 
■0(r) 



1 + 
1 - 2 



M 
r 



362 



\bQ\ 



F{a, b) is the incomplete elliptic integral of first kind. We follow the convention of and set G = 1. 
In this case one obtains the potential 



V{z) 

which can be expanded around z = as 



OO 

V{z) = 



VnZ 



(52) 

(53) 
(54) 



(55) 



(56) 
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0.01 



FIG. 6: Deflection angle for the Einstein-Born-Infeld problem as a function of ro assuming b — M = 1 and Q — 1/2. The 
potential in our formula is expanded to order z*. (color online) 



Unlike in the previous two cases V{z) is not polynomial in z and the deflection angle reads 



A(f> 



(1) 

PAIS ^ 



1 _ ML 



+ f «6 + 



(57) 



Clearly one has to keep in mind that the truncation of the series ()56|l to a finite order is an additional source of error 
in our calculation: in practice, however, it is straightforward to include further terms of the expansion. 

In Fig. we have compared the exact result obtained numerically integrating the integral in A(f> with the result 
obtained with our first order formula using the expansion V{z) to order z^: 



Viz) 



2M 
ro 



z^ + ViZ^ + 0[z^ 



(58) 



where 



7Q^ 

3rl 



I {3b' - 2/3) ^ + ( ^Q^bQ - 3bM 



y- (362-2/3) 



3 rl 



3/2 



\3VbM 



462Q4 



(59) 



and 



f3' 



7 



= F ^arccos f 1 — 



2bQ 



bQj ' V2 



(60) 
(61) 



Our analytical formula reproduces with high accuracy the numerical result obtained assuming 6 = M = 1 and 

g = 1/2. 

It is also easy to obtain the asymptotic behavior of A0 from our expression 



A0 



4M ^ ^ 24 Af 2 3t:Q'\^ ^ f 160 A-P ^ MQ 



4 



V 7r2 



O 



1 



(62) 
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V. CONCLUSIONS 



111 this paper we have presented a new method to obtain analytical expressions for the deflection angle of a ray 
of light propagating in a spherically symmetric static metric. We have been able to prove the convergence of our 
approach and to estimate the rate of convergence of the series obtained applying our method: the series converges 
exponentially and can be applied over all the physical region, as explictly seen in the case of the Schwarzschild and 
Reissner-Nordstrom metrics, where the correct location of the photon sphere is recovered. 

This method has been used to derive a first order formula, which is valid for general spherically symmetric static 
metric tensor: we have tested this formula in four different cases, observing that it is quite accurate even in proximity 
of the photon sphere. Clearly, higher order corrections to the first order formula of this paper will further improve 
the quality approximation, given the convergent nature of our expansion: we plan to study higher order corrections 
to our formula in a forthcoming paper. 

We also stress that the scries obtained with our method arc nonpcrturbativc, because they do not correspond to an 
expansion in a small parameter and therefore they arc capable of providing small errors even when the parameters in 
the model are not small (a typical perturbative parameter would be GM/vq). 
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